Statistics of Velocity Gradients in 
Two-Dimensional Navier-Stokes and Ocean Turbulence 
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Probability density functions and conditional averages of velocity gradients derived from up- 
per ocean observations are compared with results from forced simulations of the two-dimensional 
Navier-Stokes equations. Ocean data are derived from TOPEX satellite altimeter measurements. 
The simulations use rapid forcing on large scales, characteristic of surface winds. The probability 
distributions of transverse velocity derivatives from the ocean observations agree with the forced 
simulations, though they differ from unforced simulations reported elsewhere. The distribution and 
cross-correlation of velocity derivatives provide clear evidence that large coherent eddies play only 
a minor role in generating the observed statistics. 
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I. INTRODUCTION 

Statistical properties of turbulent flows, such as proba- 
bility density functions (pelfs), are important for char- 
acterizing turbulence. For instance, velocity gradients 
are directly related to velocity correlations, relative dis- 
persion, and energy dissipation in the fluid [l]]. This 
study evaluates statistics of turbulence, as observed in 
recent satellite measurements of the upper ocean. Statis- 
tics of observed phenomena are compared with corre- 
sponding statistics for the forced two-dimensional Navier- 
Stokes equations. Our results show that, in compar- 
ison with unforced decaying turbulence, simple forced 
two-dimensional Navier-Stokes equations provide better 
agreement with ocean observations. 

For this analysis, ocean velocities are derived from al- 
timeter data collected by the TOPEX/POSEIDON satel- 
lite, which performs repeated measurements of the height 
t] of the ocean surface. We use only observations from the 
TOPEX altimeter, which has lower noise levels than the 
POSEIDON instrument. The geostrophic relation, v x = 
(g/f)dr]/dy, yields the velocity component perpendicular 
to the satellite ground track. Surface geostrophic veloc- 
ities are characteristic of sub-surface flow in the ocean 
H . This geostrophic flow is typically well-represented by 
two-dimensional shallow-water equations and resembles 
two-dimensional turbulence ||, |4j. The derivative along 
the satellite track, d v v x , yields the transverse velocity 
gradient. We compute velocities v from consecutive high- 
frequency altimeter measurements jH[ ^, 0, ^), and then 
determine velocity gradients by computing along track 
differences over a distance of 12 km. For comparison, the 
first baroclinic Rossby radius ranges between 10 km and 
80 km between 60° and 10° latitude §, |§, so transverse 
gradients over 12 km distance are expected to be repre- 
sentative of mesoscale geostrophic motions. The cross- 



track, or longitudinal, derivative cannot be determined. 
Higher order derivatives are increasingly noisy. 

Earlier results have shown that velocities typically have 
Gaussian distributions within small regions of the ocean. 
When satellite data from the global ocean were com- 
bined, the resulting pdfs were non-Gaussian, due to re- 
gional variations in velocity variance @, §. When veloc- 
ities were normalized by their local variances, the pdfs 
were Gaussian B, at least for well-sampled velocities 
within three standard deviations of the mean. Similar 
results were obtained for subsurface floats deployed in 
the North Atlantic Ocean, although analysis for veloci- 
ties more than three standard deviations from the mean 



indicated non-Gaussian tails 11 1. The Lagrangian statis- 
tics of floats are however not directly comparable to the 
results from the TOPEX altimeter, which captures the 
Eulerian statistics. In this study, we specifically normal- 
ize velocities and velocity gradients by their local vari- 
ances before computing pdfs and other statistics. 

We compare observed oceanic pdfs with simulations of 
two-dimensional quasi-geostrophic flow. The equations 
of motion are 
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q = DV 2 q + F, 



(1) 



where the potential vorticity q = —^ip+ijj/R 2 . The sec- 
ond term in the potential vorticity is neglected both in 
the quasi-geostrophic limit of the shallow water equations 
and when the Rossby radius R is large in the homoge- 
neous quasi-geostrophic equations. In either case, eq. (|lj) 
is equivalent to two-dimensional Navier-Stokes flow. In 
this study, we perform simulations of the Navier-Stokes 
equations. Rapidly varying (white-in-time) forcing F is 
applied on large scales, through random stirring of the 
low vorticity modes. This forcing resembles wind forcing 
of the ocean, which varies rapidly in time but slowly in 
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space |l2], We consider an isotropic, homogeneous, 
and statistically stationary state. Simulations use a con- 
ventional pseudo-spectral method and second-order dissi- 
pation. Results were obtained on a 1024 x 1024 grid with 
long time averaging. Large-scale coherent vortices are 
clearly visible. Further details about the numerics and 
resulting velocity pdfs are described elsewhere |l4|, 15 



In many instances the velocity pdf is approximately 
Gaussian |l6), and this is also the case for the simulated 
turbulence here J^j. Far more conclusive than the ve- 
locity distribution turns out to be the statistics of veloc- 
ity derivatives. A number of authors have investigated 
the velocity gradients of three-dimensional Navier-Stokes 
turbulence || ||, [H], g || H H , but here we 
consider the far less studied two-dimensional case. Mea- 
surements of the transverse velocity derivatives are pre- 
sented in section [n| Section [II discusses evidence that 
large eddies alone provide only a minor contribution to 
the observed statistics. Section IV briefly discusses perti- 



nent differences between forced and unforced turbulence. 
The last section contains conclusions. 



II. THE PROBABILITY DISTRIBUTION OF 
VELOCITY DERIVATIVES 

Earlier work based on satellite altimeter data reported 
transverse velocity gradient pdfs in small parts of the 
ocean 0. These results differed from gradient pdfs 
derived for decaying two-dimensional turbulence, which 
showed an approximate Cauchy distribution during the 
late stage of the evolution ||(| J27|. The discrepancy is 
resolved by comparing to simulations of stationary tur- 
bulence. 

Figure [l] shows velocity gradient pdfs derived from 
ocean observations and simulations. The solid line in- 
dicates the pdf of normalized velocity gradient data de- 
rived from global satellite altimetry. To determine the 
oceanic pdf, velocity gradient data drawn from latitudes 
between 10° and 60°N and between 10° and 60°S are 
sorted geographically into 2.5° by 2.5° boxes. Data near 
the equator are omitted because the geostrophic relation- 
ship is weak at low latitudes. The standard deviation 
of gradients in each latitude-longitude box varies from 
1.4 x 10" 5 s _1 near 60° latitude up to 6.4 x 10 -5 s _1 near 
10° latitude, with a median value of 1.9 x 10~ 5 s _1 . To 
compensate for this geographic variation, gradients are 
normalized to have unit standard deviation in each box, 
and then the pdf is computed from all of the normalized 
gradient data. For comparison, we also normalized our 
pdfs using the mean absolute value of the velocity gradi- 
ent; this did not diminish the strong tails of the gradient 
pdf. 

The dashed line in Fig. [t] represents the transverse ve- 
locity gradient pdf from two-dimensional Navier-Stokes 
turbulence. The dotted lines represent the narrow Gaus- 
sian distribution and the broader Cauchy distribution, 
P(x) = c/[tt(c 2 + x 2 )}, with long tails. The tails con- 
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FIG. 1: (a) Global variance-normalized pdf of the veloc- 
ity gradient in the ocean (solid line) compared with simu- 
lations of two-dimensional Navier-Stokes turbulence (dashed 
line) on a semi-logarithmic scale, (b) Same quantities on a 
linear scale. In both Gaussian and a Cauchy distri- 

bution are shown for comparison (dotted lines). The ocean 
pdf is averaged over about 13 million data points. Data are 
normalized by the standard deviation a, as described in the 
text. 



tribute noticeably to the standard deviation of the pdf, 
and therefore the Gaussian is fitted to the data without 
requiring unit area and unit standard deviation. This is 
necessary to make the Gaussian closely approximate the 
central part of the pdf. Since the Cauchy distribution 
cannot be normalized by its standard deviation, the con- 
stant c is chosen such that P(0) matches. Both the sim- 
ulated and observed gradient pdfs appear Gaussian for 
small velocity gradients, up to about one standard devi- 
ation. For large gradients they decay significantly more 
slowly than do Gaussian tails but substantially faster 
than the Cauchy distribution found in simulations of de- 
caying turbulence |7). There is good agreement be- 
tween observed and simulated pdfs up to even the largest 
fluctuations measured in the simulation. 

Error bars for the pdfs were estimated by grouping the 
data into N groups and computing pdfs for each group. 
Error of the mean pdf is then taken to be the standard 
deviation of the pdf divided by y/N. For this analysis, N 
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was the total number of 2.5° boxes for the surveyed ocean 
or the number of velocity snapshots for the simulation. 
Since many ocean observations are available, statistical 
errors are expected to be small compared to systematic 
errors. In fact, the statistical errors are frequently nar- 
rower than the line width in Fig. |l|. Differences between 
the two distributions exceed the statistical errors and are 
likely to be due to a number of factors, including instru- 
mental and atmospheric correction errors in the altimeter 
data, which make the measurements noisy, as well as dif- 
ferences in the physics of two-dimensional Navier-Stokes 
equations compared with the ocean. 

The kurtosis (flatness), (x 4 ) / (x 2 ) 2 , can serve as quan- 
titative comparison of the shape of the pdf. In the sim- 
ulation results, the velocity gradient pdf has a kurtosis 
of 4.7, indicating clear deviation from Gaussian distribu- 
tion. If the observed ocean pdf is terminated beyond the 
extent of the simulated one, at about four and a half stan- 
dard deviations, its kurtosis is also 4.7. This quantitative 
comparison confirms that the oceanic pdf is substantially 
better matched by the simulation than by either of the 
two ideal distributions. 

Velocity gradient pdfs depend on the spatial separa- 
tion between velocity measurements. The velocity corre- 
lations between two points decrease with distance, and 
velocities at points very far apart can be assumed to be 
statistically independent. The distribution of velocity 
differences across very large distances reduces to that of 
the velocity (with twice the variance). The 12 km separa- 
tion of TOPEX observations is small compared with the 
decorrelation length scales of wind forcing, 0(1000 km), 
and of mesoscale ocean features 0(100 km), so veloci- 
ties at adjacent observation points are expected to be 
strongly correlated. Therefore, to obtain comparable re- 
sults from the numerical simulation, we have computed 
gradient pdfs from velocity differences over asymptoti- 
cally small separations. For comparison, if we compute 
gradient pdfs over progressively larger distances in the 
simulation, then the distribution narrows from its orig- 
inal shape (dashed line in Fig. |l|) and becomes close to 
Gaussian. 

The basic simulations had a large-scale Reynolds num- 
ber on the order of 10 , while for ocean turbulence a 
Reynolds number of 10 7 might be typical |2q] . Pdfs were 
also determined for simulations with lower and higher 
Reynolds numbers, using respectively lower and higher 
resolutions, but shorter sampling time. There is no signif- 
icant change in the shape of the pdfs 1 29 , although these 
data do not exclude a weak dependence on Reynolds 
number. The absence of any detectable Reynolds number 
dependence suggests that the simulation data are close 
to what they look like at substantially higher Reynolds 
number. The difference in the length of the tails in Fig. |l| 
may be due to the vast difference in Reynolds number, 
difference in sampling size, and errors from the numerical 
differentiation of data. 

The real ocean differs from the forced Navier-Stokes 
system because of the addition of the (3-eSect, stratifica- 
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FIG. 2: Pdf of the velocity gradient in simulations of two- 
dimensional Navier-Stokes turbulence. The solid line is for 
the velocity gradients produced by large coherent vortices. 
The dashed line corresponds to the complete flow field. Both 
distributions are normalized by the same standard deviation, 
hence preserving differences in their width. 



tion, three-dimensional motions, and buoyancy. Hence it 
is surprising that there is such a close agreement between 
measurement and simulation. In any case, the agree- 
ment between observation and simulation suggests that 
the oceanic velocity statistics may be understood within 
the framework of two-dimensional turbulence. 



III. THE ROLE OF COHERENT VORTICES 

Idealized models of point vortices predict a Cauchy dis- 
tribution for the velocity gradients and a Gaussian dis- 
tribution for the velocity |^6|, ^7], [$0| |3l], |3^] . This agrees 
with results from decaying two-dimensional turbulence 
p6j Hfl . Hence, in the late stages of decay, the statis- 
tics of velocity gradients have been successfully under- 
stood to result from the far-field of well-separated vor- 
tices 1 26, p7|. In contrast, pdfs of ocean surface veloc- 



ity gradients are observed to have more rapidly decaying 
tails than do Cauchy distributions. Also as Fig. [I] makes 
evident, the simulations of stationary two-dimensional 
turbulence show far less pronounced tails than a Cauchy 
distribution. 

The discrepancy arises not only in the shape of the 
distribution but also in its width. A straight-forward 
way of illustrating this is to calculate the velocity field 
produced by vortices with vorticities that exceed twice 
the root-mean-square value. Figure ^| shows the trans- 
verse velocity gradients produced by these coherent vor- 
tices (solid line) . For comparison, the actual distribution 
is shown as a dashed line. Clearly the large coherent 
vortices do not generate enough intermediate gradients. 
(Nor, for that matter, do they account for most of the 
velocities.) Consequently, the distribution of gradients is 
poorly accounted for by large-scale coherent vortices. 

The contribution of the small-scale turbulence is also 
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FIG. 3: Conditional average of transverse velocity gradient 
with velocity, yj {(d y v x ) 2 \ \v x \), for two-dimensional Navier- 
Stokes (dashed line with error bars) and ocean turbulence 
(solid line). The conditional average produced by the large 
vortices in the simulation is also shown (dash-dotted line). 
All three graphs are normalized by their respective standard 
deviations. The thick solid line indicates growth proportional 
to the velocity. The error bars include only the statistical 
error expected from averaging of the 32 snapshots, showing 
twice the standard error of the mean. The inset shows the 
conditional average of the absolute value { | c?^ | | \v x \). 



relevant. This agrees with the basic physical picture, 
according to which the late stage of decaying turbulence 
consists of coherent vortices. Its statistics can therefore 
be understood in terms of them. In the stationary case, 
on the other hand, fluctuations over a wide spectrum of 
spatial scales contribute to the gradient statistics. 

Available statistical variables from the altimeter are 
the velocity and the transverse velocity derivative. 
Hence, one can study the cross-correlation between these 
two quantities. Here, we examine the conditional av- 
erage of the squared velocity gradient as a function of 
velocity, ({d y v x ) 2 \ \v x \), which is the average of (d y v x ) 2 
over all points with velocity component ±v x . The slope 
of ((d y v x ) 2 \ \v x \) is a measure of the correlation between 
the velocity and the transverse velocity derivative. If 
there were no correlation between the velocity at a point 
and the gradient at the same point, the conditional av- 
erage would be constant for all values of v x and would 
be exactly one if velocity gradients were normalized by 
their standard deviation. In contrast, if gradients and 
velocities were strongly correlated, as would be expected 
around an isolated vortex, then the graph for the condi- 
tional average would have a pronounced slope. 

Figure || shows the square-root of the measured con- 
ditional average together with that for two-dimensional 
Navier-Stokes turbulence. For the ocean, we have nor- 
malized both v x and d y v x by their standard deviations 



in each 2.5° by 2.5° geographic box, because both quan- 
tities vary spatially. The axes are labeled in units of their 
respective standard deviations, \J {v 2 } and \J ((d y v x ) 2 ). 
The correlation between velocity and its gradient is weak, 
but the gradients tend to be higher when the velocity is 
large. If oceanic gradients beyond four standard devia- 
tions are excluded, which is a fairer comparison with the 
simulation, the conditional average is closer to one. The 
longitudinal component of the conditional average (not 
shown) exhibits behavior similar to the transverse com- 
ponent. Also shown in Fig. || is the conditional average 
for the velocity field of vortices larger than twice the root- 
mean-square vorticity (dash-dotted line). As expected 
there is a comparatively strong correlation between ve- 
locities and velocity derivatives. At large velocities the 
slopes of the graphs for the coherent vortices and the 
ocean are similar. This may indicate influence by large 
eddies in regions where the velocities are high, although 
no corresponding evidence is found in the probability dis- 
tribution of the gradients. The situation at high veloci- 
ties is therefore somewhat ambiguous. For small veloci- 
ties, which cover most of the area, conditional averages 
are near one for simulations (dashed line) and observa- 
tions (solid line), indicating that at low velocity, gradi- 
ents arc almost uncorrelated with velocity. The inset in 
Fig. H shows the conditional average using Idyi^l instead 
of (d y v x ) 2 , which is less sensitive to outliers. For small ve- 
locities, the agreement between observation and simula- 
tion is closer and the discrepancy between the large eddy 
field and the other two conditional averages is stronger. 

The deviation of the observed conditional averages 
from that for coherent vortices strengthens the evidence 
that the gradient statistics are unaccounted for by the 
velocity field created by large-scale coherent eddies. The 
role of coherent vortices in generating the observed veloc- 
ity statistics is minor. This conclusion cautions against 
attempts to model oceanic velocity fields by large eddies. 



IV. FORCED VERSUS UNFORCED 
TURBULENCE 

Although only the transverse velocity component can 
be determined from altimeter data, simulations also per- 
mit us to examine the longitudinal derivative, d x v x . 
Figure |^ shows a clear difference between the behav- 
ior of the longitudinal and transverse components. In 
our forced simulations, the standard deviation of lon- 
gitudinal fluctuations is about 60% of that for trans- 
verse fluctuations. In isotropic and incompressible tur- 
bulence there is an exact relation between the standard 
deviation of transverse and longitudinal component [ p4[ . 
With a calculation analogous to the well-known three- 
dimensional case p4[ , we find in the two-dimensional case 
{{d y v x ) 2 ) = 3 {(d x v x ) 2 ). Hence, the standard deviation 
for the longitudinal component is l/-\/3 ~ 58% of that 
for the transverse component. This agrees with the mea- 
sured value of 60% in the simulation. 
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(b) d x v x> d y v x 

FIG. 4: (a) Probability density functions of velocity deriva- 
tives for forced two-dimensional Navier-Stokes turbulence on 
a semi-logarithmic scale, (b) Same quantities on a linear scale. 
In both panels, the dashed line shows the transverse compo- 
nent, dyVx, and the solid line the longitudinal component, 
d x v x . Both are normalized by the standard deviation of the 
transverse component a — -J ((d y v x ) 2 ). The dotted line is 
a Gaussian. Transverse and longitudinal gradient pdf differ 
from each other in width and shape. 



frequently than do large transverse gradients. For sim- 
ple point-vortex models both, transverse and longitudi- 
nal components, are distributed like Cauchy distributions 
(albeit with different standard deviations) [E(J. Also in 
the late stage of decaying turbulence, the transverse com- 
ponent is distributed in the same way as the longitudinal 
component |2(|. This is yet another difference between 
forced and unforced turbulence. 

Overall, our study establishes a clear distinction be- 
tween the gradient statistics of unforced (freely decay- 
ing) and forced (stationary) turbulence. The presence of 
forcing not only influences the properties of large-scale 
vortices but also changes the distribution of eddies over 
different scales. (Freely decaying turbulence has an in- 
verse energy cascade while two-dimensional turbulence 
forced at large scales is governed by a direct enstrophy 
cascade.) Further study is needed to determine how the 
statistics may depend on the temporal and spatial struc- 
ture of the forcing. 



V. CONCLUSIONS 

In conclusion, we find that transverse velocity deriva- 
tive pdfs from observed upper-ocean turbulence agree 
closely with forced two-dimensional simulations but dif- 
fer from previously reported unforced turbulence. The 
forcing diminishes the role of coherent vortices in the per- 
tinent statistics. The distribution and cross-correlation 
of velocity derivatives provide clear evidence that large 
coherent eddies play only a minor role in generating 
the observed statistics. Further study of forced two- 
dimensional turbulence appears likely to shed light on 
the character of meso-scalc turbulence in the ocean. 
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